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Abstract. We present a new determination of the metallicity distribution, age, and luminosity function of the 
Galactic bulge stellar population. By combining near-IR data from the 2MASS survey, from the SOFI imager 
at ESO NTT and the NICMOS camera on board HST we were able to construct color-magnitude diagrams 
(CMD) and luminosity functions (LF) with large statistics and small photometric errors from the Asymptotic 
Giant Branch (AGB) and Red Giant Branch (RGB) tip down to ~ 0.15Mq. This is the most extended and 
complete LF so far obtained for the galactic bulge. Similar near-IR data for a disk control field were used to 
decontaminate the bulge CMDs from foreground disk stars, and hence to set stronger constraint on the bulge age, 
which we found to be as large as that of Galactic globular clusters, or >10 Gyr. No trace is found for any younger 
stellar population. Synthetic CMDs have been constructed to simulate the effect of photometric errors, blending, 
differential reddening, metallicity dispersion and depth effect in the comparison with the observational data. By 
combining the near-IR data with optical ones, from the Wide Field Imager at the ESO/MPG 2.2m telescope, a 
disk-decontaminated (Mk, V — K) CMD has been constructed and used to derive the bulge metallicity distribution, 
by comparison with empirical RGB templates. The bulge metallicity is found to peak at near solar value, with 
a sharp cutoff just above solar, and a tail towards lower metallicity that does not appreciably extend below 
[M/H]~ -1.5. 
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1. Introduction 

Galactic spheroids, i.e., elliptical galaxies and the bulges 
of spiral galaxies, contain a large fraction, probably the 
majority, of all the stellar mass in the local universe 
(Persic & Salucci 1992; Fukugita, Hogan, & Peebles 1998). 
Understanding their formation and evolution is therefore 

* Based on observations collected at the European Southern 
Observatory, La Silla, Chile, obtained from the ESO/ST- 
ECF Science Archive Facility, and on observations with the 
NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, operated by AURA Inc. under con- 
tract to NASA 



crucial to understand galaxy formation in general, and 
to reconstruct the whole cosmic star formation history. 
However, it is only in the bulge of our own Galaxy that 
we can resolve the stellar population all the way to the 
bottom of the main sequence (MS), construct accurate 
color-magnitude diagrams (CMD) and luminosity func- 
tions (LF), and then address a series of issues of great im- 
portance for a better understanding of galactic spheroids 
in general. Such issues include the direct measure of some 
of the fundamental ingredients of galaxy models, like the 
initial mass function (IMF), the mass-to- light ratio and 



2 



M. Zoccali et al.: Age and Metallicity Distribution of the Galactic bulge 



quantities such as the stellar age and metallicity distribu- 
tions (MD). 

Efforts in these directions have been quite numerous in 
recent years, and reviewing all of them is well beyond the 
scope of this paper. However, here are briefly mentioned 
some of the most recent attempts at determining the age, 
IMF and MD of the bulge stellar population. 

The only published study of the bulge MD based 
on high resolution (R=17,000) spectra remains that by 
McWilliam & Rich (1994; hereafter M&R). They analyzed 
a sample of 11 red giants in Baade's Window, and used 
these data to re-calibrate a previous bulge MD based on 
low resolution spectra for 88 stars (Rich 1988). The result- 
ing MD is centered around [Fe/H]=— 0.2, with some 34% 
of the stars above solar metallicity, and no stars below 
[Fe/H]~ -1.3. 

The bulge MD has also been derived using lower reso- 
lution spectra. Sadler et al. (1996) measured line strength 
< Fe > indices in low resolution (R ~ 1000) spectra of 
a sample of 322 bulge giants. These indices were cali- 
brated vs. [Fe/H] following Faber et al. (1985). Although 
the derived [Fe/H] abundances for the few stars in com- 
mon with M&R are in good agreement, the global MD by 
Sadler et al. (1996) is symmetric around [Fe/H]=0, and 
shows a population of stars with metallicity above solar 
more prominent than in any other study, which may re- 
sult from having applied the <Fe> vs. [Fe/H] calibration 
beyond its range of validity. More recently, Ramirez et al. 
(2000) derived the bulge MD by measuring the equivalent 
width of Ca, Na and CO lines in a sample of 110 M giants 
observed with resolutions from R^ 1300 to 4800. They 
converted equivalent widths into [Fe/H] by means of the 
Frogel et al. (2001) calibration, based on Galactic glob- 
ular clusters (GCs) of independently known [Fe/H]. This 
assumes that the a-element enhancement of the bulge is 
the same as that of the clusters. The MD obtained in this 
case is similar to that of M&R, but with a sharper peak 
at [Fe/H]=0. The appreciable differences among the MDs 
derived in these studies are likely to arise from uncertain- 
ties in the calibration of the low resolution indices, that 
may introduce important systematics. 

Dating bulge stars is complicated by several factors, 
such as crowding, depth effects, variable reddening, metal- 
licity dispersion, and contamination by foreground disk 
stars. From WFPC1 observations of the bulge field known 
as Baade's Window, Holzman et al. (1993) inferred a dom- 
inant intermediate age component in the Galactic bulge. 
Shortly after the refurbishment of HST, two metal rich 
globular clusters of the bulge (NGC 6528 and NGC 6553) 
were observed with WFPC2 (Ortolani et al. 1995, here- 
after Paper I). These clusters are respectively at ~ 4° and 
~ 6° from the galactic center, and their overall metal- 
licity [M/H] is about solar (Barbuy et al. 1999; Cohen 
et al. 1999), close to the average for stars in Baade's 
Window (M&R). Like most other clusters within ~ 3 kpc 
from the Galactic center, they belong to the population of 
bulge globular clusters that have age, kinematical proper- 
ties and metallicity distribution that are indistinguishable 



from those of bulge stars (e.g., Minniti 1995; Paper I; Cote 
1999). Hence, we do not hesitate to refer to them as bulge 
globular clusters (see also Forbes, Brodie & Larsen 2001; 
Harris 1976, 2001; Zinn 1996). In the case of NGC 6528 
and NGC 6553 the bulge membership is also confirmed 
by the proper motion of the two clusters, in both cases 
well within the proper motion distribution of bulge stars 
(Zoccali et al. 2001a; Feltzing, Johnson, & de Cordova 
2002). In Paper I it was shown that i) the two clusters have 
virtually identical CMDs, indicating that they have the 
same age and metallicity; ii) their magnitude difference 
between the horizontal branch (HB) clump and the MS 
turnoff is virtually identical to that of the inner halo glob- 
ular cluster NGC 104 ([Fe/H]=-0.7), and in) the LF of 
NGC 6528 (the least reddened of the two clusters) is virtu- 
ally identical to the LF of bulge stars in Baade's Window, 
when allowance is made for the distribution of bulge star 
distance along the line of sight. From all these evidences 
Ortolani et al. concluded that i) the bulge underwent rapid 
chemical enrichment to solar abundance and beyond, very 
early in the evolution of the Milky Way (MW) Galaxy; ii) 
the bulk of bulge stars formed nearly at the same time as 
the halo globular clusters, and Hi) no more than ~ 10% 
by number of the bulge population can be represented by 
intermediate age stars. These conclusions have been fur- 
ther strengthened after a statistical decontamination of 
the bulge CMD from the foreground disk stars (Feltzing 
& Gilmore 2000). 

Dating the bulge via the MS turnoff also allows to 
check the usefulness of asymptotic giant branch (AGB) 
stars brighter than the red giant branch (RGB) tip as age 
indicators. Indeed, being such AGB stars much brighter 
than the MS turnoff, they have often been used in the 
attempt to date other galactic spheroids (such as M32 
and the bulge of M31), sometimes inferring intermediate 
ages for some of these systems (e.g, Davidge and Nieto, 
1992; Freedman 1992; Elston and Silva, 1992). However, 
the existence of AGB stars brighter than the RGB tip 
does not ensure a population to be of intermediate age: 
the old, metal rich ([Fe/H]> — 0.7) globular clusters of 
the MW bulge and inner halo (such as, e.g., NGC 6553 
and NGC 104) do indeed contain some AGB stars that 
are ~ 1 magnitude brighter than the RGB tip (Frogel 
and Elias 1988; Guarnieri, Renzini, & Ortolani 1997). In 
any event, the RGB+AGB luminosity function conveys 
important information concerning the stellar population 
of the bulge, and its better understanding is important for 
the interpretation of the observations of nearby spheroids. 

At the opposite extreme of the luminosity range, the 
LF of the lower MS is the only part of the LF that de- 
pends on the IMF, and therefore it allows the determi- 
nation of this important property of stellar populations. 
Recent near-infrared data obtained with HST+NICMOS, 
have provided the deepest ever J, H photometry in the 
bulge, and allowed the determination of the IMF from 
~ 1M Q down to ~ 0.15M© (Zoccali et al. 2000, hereafter 
Paper II). The IMF resulted to be quite flat, with a slope 
-1.33 ± 0.07 (compared to -2.35 for the Salpeter IMF), 
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with a hint of a steepening toward the high mass end (see 
also Holzman ct al. 1998). 

The study of the NICMOS field extended to only 408 
square arcsec [], and therefore the number of stars on the 
upper MS is too small to allow any reliable location of 
the MS turnoff, while the evolutionary stages beyond the 
turnoff are not sampled at all. Thus, in Paper II an at- 
tempt was made to extend the main sequence LF to the 
evolved stars, using near-infrared data from the litera- 
ture. A complete LF, extending from the bottom of the 
MS all the way to the AGB, was constructed by match- 
ing the NICMOS luminosity function of the lower MS, to 
the Tiede, Frogel & Terndrup (1995) LF extending from 
slightly above the turnoff to slightly above the HB, and 
finally to the Frogel & Whitford (1987) LF for the bulge 
M giants above the HB. While the resulting composite 
LF was constructed with the best data available in each 
luminosity range, several limitations were also obvious. 
Specifically, the MS turnoff region was not well sampled 
by neither the NICMOS data from Paper II (too few stars 
in the small field), nor by the Tiede et al. LF (not deep 
enough to reach the turnoff). Moreover, the upper RGB 
and AGB from Frogel & Whitford (1987) included only 
the M stars selected from an objective prism survey, hence 
omitting any luminous K giants that may be present in 
the explored bulge field. Finally, the NICMOS observa- 
tions were conducted in a field 6° from the Galactic cen- 
ter, while the two other datasets were relative to Baade's 
Window, at 4° from the center. 

In this paper we present a thorough attempt at over- 
coming the current limitations of the available CMDs and 
LFs of the bulge, thereby producing state of the art CMDs 
and LFs used for new determinations of the age and metal- 
licity distribution of the bulge. This approach is based 
on new near-IR (J, H, K s ) data obtained with SOFI at 
the ESO New Technology Telescope (NTT), as well as on 
optical (V, I) data obtained with the Wide Field Imager 
(WFI) at the ESO/MPG 2.2m telescope. These data are 
then coupled with the NICMOS data from Paper II for 
the lower MS, and with the 2MASS near-IR data for the 
upper RGB and AGB (van Dyk, 2000; Carpenter, 2001; 
and references therein). 

The paper is organized as follows: Observations and 
data reduction procedures are presented in Section 2, the 
CMDs are displayed and discussed in Section 3, while 
Section 4 is devoted to the derivation of the metallicity 
distribution. The luminosity function of the bulge is con- 
structed and discussed in Section 5, while Section 6 is de- 
voted to determining the age of the bulge, and Section 7 to 
a discussion of the RGB tip in the different bands. Finally, 
the main results are again summarized in Section 8, along 
with some inferences on the formation of the Galactic 
bulge and spheroid, also in the context of the current ev- 
idence on the formation of galactic spheroids in general. 



1 The size of the NICMOS field was erroneously quoted as 
506 square arcsec in Paper II. 



2. Observations and Data Reduction 

The data presented here come from the combina- 
tion of near-IR J, H, K s observations conducted with 
SOFIQNTT, implemented with public data from the 
2MASS survey and with NICMOS data from Paper II, 
plus optical V,I images taken with the WFI@2.2m as 
part of the EIS PRE-FLAMES survey (Momany et al. 
2001; and references therein), and retrieved from the ESO 
archive. The relevant technical information about the ob- 
servations are reported in Table [j], while details on the 
observational strategy and data reductions are described 
here below. 

2.1. Near-IR data 

A bulge mosaic field centered at (1, 6) = (0.277, —6.167) 
(RA=18:11:13, DEC=-31:43:49; J2000) was observed 
with SOFI@NTT, through the filters J,H,K S , in the 
nights of 9-12 June 2001. This particular field was selected 
for the present study because it includes the field already 
covered by the NICMOS observations (Paper II), which in 
turn was chosen for having a reddening (E{B — V) = 0.47) 
as low as that of the most widely studied Baade's Window 
while being less crowded. Indeed, being ~ 2° further away 
from the Galactic center the density of stars is lower by a 
factor ~ 2, which allows more accurate photometry espe- 
cially at very faint magnitudes. 

In order to optimize the photometry of both the turnoff 
region and the upper RGB, the observations were split 
in short and long exposures, each with a different field 
objective, and therefore mapping areas of different sizes. 
For a good sampling of RGB+AGB stars, short exposures 
were obtained for a 8'3 x 8.'3 region, with a mosaic of 
four fields of the low resolution camera, which has a pixel 
size of 0'.'292 and field of view of 4f9 x 4.'9 (hereafter SOFI- 
LARGE field). In order to better sample the PSF, a second 
mosaic of four deeper exposures mapping a smaller area 
(3.'9 x 3'9) was secured with the same camera coupled with 
the focal elongator, therefore yielding a pixel size of 0V144 
and a 215x215 field of view (hereafter SOFI-SMALL field). 
An overlap of ~ 20% was always present among the four 
fields of each mosaic. 

Exposures of a disk control field located at (I, 6) = (30,0) 
(RA=19:07:32, DEC=-05:19:57) were also obtained for 
one pointing of the large field camera in order to estimate 
the disk contribution to the bulge CMD (Section 0). 

Each of the deep exposures, both for the disk and 
the bulge field, was obtained with detector setup DIT=6 
NDIT=5 (i.e, every frame was the average of 5 exposures 
of 6 sec. each), while for the shallow ones we used DIT=1.2 
and NDIT=5^|. Both of them were repeated with a ran- 
dom dithering pattern until reaching the exposure time 
listed in Table [|. Both sky transparency and seeing were 
somewhat variable during the run, and for this reason the 

2 See the SOFI User Manual for a more detailed explanation 

of the instrument setups: 

|http : //www . eso . org/las ilia/Telescopes /NEWNTT/s of i 
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Table 1. Log of the observations 



Date 


Object 


Camera 


RA 


DEC 


Filter 


Exptimc 


11 Jun 2001 


bulge 


SOFI-LARGE 


18:11:05.0 


-31:45:49 


.1 


36s 


// 


// 


// 


II 


a 


H 


36s 












Ks 


36s 








18:11:22.0 


-31:45:49 


.1 


36s 


// 


// 


a 


II 


II 


H 


36s 


// 


a 


a 


II 


II 


Ks 


36s 


// 




a 


18:11:22.0 


-31:41:49 


.1 


36s 


// 


a 


a 


II 


II 


H 


36s 


// 


a 


a 






Ks 


36s 


// 




a 


18:11:05.0 


-31:41:49 


.1 


36s 


// 


a 


a 


II 


// 


H 


36s 


// 


a 


a 


II 


II 


Ks 


36s 


11-12 Jun 2001 


bulge 


SOFI-SMALL 


18:11:05.0 


-31:45:49 


.1 


1180s 




// 


// 


II 




H 


1530s 












Ks 


2880s 


9-11 Jun 2001 




// 


18:11:13.5 


—31:45:49 


.1 


1200s 




// 


II 




II 


H 


1020s 












Ks 


1920s 


9-11 Jun 2001 




If 


18:11:13.5 


-31:43:49 


.1 


1200s 




// 


II 






H 


1020s 






II 






Ks 


1920s 


9-11 Jun 2001 




If 


18:11:05.0 


-31:43:49 


.1 


1200s 


// 


// 


II 


II 


II 


H 


1020s 


// 


// 


II 


II 


II 


Ks 


1920s 


11 Jun 2001 


disk 


SOFI-LARGE 


19:07:32.0 


-05:19:57 


J 


600s 


// 


// 


// 


// 


// 


H 


510s 


// 


// 


II 


// 


// 


Ks 


960s 


15 Apr 1999 


bulge 


WFI 


18:10:17.0 


-31:45:16 


V 


20s 


// 


// 


WFI 


// 


// 


V 


2 x 300s 


// 


// 


WFI 


// 


// 


I 


20s 


// 


// 


WFI 


II 


// 


I 


3 x 300s 



total exposure times on different fields were adjusted in 
order to compensate for these effects. The typical seeing 
during the observations was f//80 ± 0'.'15. 

Images were pre-reduced using standard IRAF rou- 
tines. A sky image, obtained by median combination of 
the dithered frames of each filter, was subtracted from 
each frame. Flat fielding was then performed using the 
" SpecialDomeFlat" template which applies the appropri- 
ate illumination corrections, as described in the SOFI User 
Manual. Finally, all the dithered frames obtained in a se- 
quence were averaged in a single image for each filter. In 
what follows we will call "frame" the combination of each 
of these sets. 

Standard photometry, including PSF modeling, was 
carried out on each frame using the DAOPHOTII pho- 
tometry package (Stetson 1987). We used all the stars 
identified in each frame to obtain the coordinate transfor- 
mations among the frames. These transformations were 
used to register the frames and obtain a median image. 
The latter, having the highest S/N, was used to create 
the star list, by means of a complete run of DAOPHOTII 
and ALLSTAR. The final star list, together with the co- 
ordinate transformations, was finally used as input for 
ALLFRAME (Stetson 1994), for the simultaneous PSF- 
fitting photometry of all the frames of each field. 

Only one of the four half-nights assigned to this pro- 
gramme was photometric. During that night we observed 
several standards from the Persson et al. (1998) catalog, 



three of which just before and after one of the bulge fields 
observed with SOFI-SMALL. This field was calibrated 
by means of the standard stars, whose frames were pre- 
processed in the same way as the science ones. Aperture 
photometry within a radius of 5.8 arcsec was obtained for 
the standard stars, and aperture corrections determined 
using some bulge isolated and unsaturated stars were ap- 
plied to the PSF instrumental magnitudes of the bulge 
field. The calibration equations obtained from the stan- 
dard stars were then applied to the magnitude of the bulge 
stars, neglecting the color term due to the very small color 
range of the near-IR standard stars. Having calibrated one 
bulge field, all the others were registered to the same pho- 
tometric system by comparison of the common stars: a 
single large field included the whole mosaic of four small 
fields, and there was always ~ 20% overlap among the 
four fields of each mosaic. 

The whole bulge area mapped with SOFI is also in- 
cluded in the 2MASS near-IR survey (Carpenter 2001), 
whose second incremental release is publicly available on 
the WEB0. Figure |l] shows the comparison between the 
J, H , K s magnitudes from our calibrated photometry and 
the same quantities for the stars in common with the 
2MASS point source catalog. The latter has a limit mag- 
nitude of K s ~ 15, therefore only the stars measured in 
one of the four SOFI-LARGE fields (i.e., in the short ex- 



tittp : //irsa. ipac . caltech. edu/applications/2MAS£ 
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Fig. 1. Comparison between the SOFI photometry and 
that from the 2MASS survey, for the stars in common. 
The solid line is the least square fit to the data, while the 
dotted one is the relation found by Carpenter (2001). 



Fig. 2. Comparison between the SOFI photometry pre- 
sented here and the deep NICMOS photometry from 
Paper II, for the stars in common. Open circles refer to 
stars that may be affected by blending in SOFI frames. 



posures) are plotted in this figure. The solid line repre- 
sents a least square fit to the data, while the dotted one 
is the relation between the 2MASS and the Persson et al. 
(1998) photometric system (Carpenter 2001). The latter 
is based on 2MASS observations of 82 standard stars from 
Persson's catalog. 

A trend with magnitude is clearly visible in Fig. [l], 
with the difference between the two photometries being 
zero at the faintest limit but increasing up to ~ 0.1 mag 
at the brightest end in J and H , while remaining quite 
small (<0.05 mag) in K s . Such a behavior points to a 
non-linearity effect in one of the two detectors. To our 
knowledge neither the 2MASS nor the SOFI measure- 
ments should be affected by non-linearity; in particular 
the SOFI detector has been tested to be linear within 2% 
up to 10,000 ADU (cf. the SOFI User Manual), and only 
the three brightest stars in this plot have counts above this 
limit. We are then left with no explanation for the trend 
seen m Fig. H but none of our conclusions relies on an 
accuracy better than ~ 0.1 in the magnitudes. Moreover, 
if the error is in our measurements and not in the 2MASS 
ones, the trend is identical in the J and H bands, therefore 
it would affect the magnitudes but not the J — H color we 
used to construct the CMD shown in Section ||. 

The disk control field was not observed during the pho- 
tometric night, therefore the only way we had to calibrate 
its CMD was to rely on the comparison with 2MASS. 
However, given that our main interest was to put the 
disk stars in the same photometric system as the bulge 
ones, we determined the calibration transformations be- 



tween our disk instrumental magnitudes and 2MASS, but 
then also applied the small differences shown in Fig. |l|, in 
order to try and keep consistency between the disk and 
bulge calibrated photometry. 

One of the goals of this project is to obtain a com- 
plete luminosity function for the galactic bulge, from the 
RGB tip down to the faintest stars (0.15M©) measured 
with NICMOS. Therefore, it is crucial to make sure that 
the NICMOS photometry is consistent with that obtained 
here with SOFI, especially given the large uncertainties in 
the NICMOS calibration (Stephens et al. 2000). Figure | 
shows the differences in J and H for the stars in common 
between the two data sets. Given the large difference in 
spatial resolution between SOFI and NICMOS, this com- 
parison may be affected by systematics due to the pres- 
ence of close pairs of stars resolved by NICMOS but not 
by SOFI. For this reason, different symbols (filled) were 
adopted in Figure ^| for the stars whose positions matched 
within 0.5 SOFI pixels. These stars are less likely to be 
blends, because the presence of a companion resolved by 
NICMOS would likely result in a displacement of the star 
centroid. The median differences, in J and H, between the 
common stars, if only the bona fide single stars are con- 
sidered, are AJ = —0.006 and AH = —0.04 as shown in 
the figure labels. If all the common stars were considered, 
the median differences would be anyway negligible for our 
purposes: AJ = AH = -0.09. 

Completeness estimates were determined via artificial- 
star experiments. A total of about 3500 stars were added 
to the original SOFI frames, both for the bulge and for the 
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Fig. 3. Difference between the input and the output mag- 
nitude of the artificial star experiments on the SOFI 
frames. The results obtained for the SOFI-LARGE and 
SOFI-SMALL frames were matched at J = 16 as in the 
photometry of the real stars. 



disk fields, with magnitudes and colors consistent with the 
RGB+MS instrumental fiducial lines. In order to avoid to 
artificially increase the crowding, at the same time opti- 
mizing the CPU time, in each independent experiment the 
artificial stars were added along the corners of an hexag- 
onal grid, as explained in Paper II. As usual, photometry 
of the artificial frames was performed following the same 
procedures as for the original ones. 

Figure || shows the results of these experiments, as the 
difference between the input and output magnitude of the 
artificial stars in each filter. 

Note that the distribution is asymmetric about the 
zero, as occasional blendings result in brighter output 
magnitudes compared to the input ones. However, the 
ridge line (i.e., the peak of the distribution) remains close 
to zero up to J ~ 18, and therefore the estimate of the 
turnoff magnitude, from the CMD ridge line is not affected 
by blending. We expect instead the luminosity function to 
be affected (see Section 5.1). The simulations have shown 
that the SOFI photometry is more than 50% complete 
above J ~ 19 and H = K* ~ 18. 



2.2. Optical data 

A 34 x 33 arcmin bulge field including the whole area 
mapped with SOFI was observed with the wide field 
imager WFI@2.2m telescope, as part of the EIS PRE- 
FLAMES programme (Momany et al. 2001). The V and 
/ raw frames have been retrieved from the public ESO 



archive, while the pre-processed images were not released 
yet. The frames were taken under very good seeing con- 
ditions, with a PSF FWHM of 0"6, as measured on 
the frames. Details about the observations are given in 
Table [l]. Science images were de-biased and flatficlded by 
means of standard IRAF routines, using a set of sky flat- 
fields taken the same night. 

The photometry for each of the 8 WFI chips was per- 
formed separately, in order to properly model the PSF 
variation across them. As usual, complete photometry and 
PSF modeling were carried out on each of the I and V 
frames, and then coordinate transformations were calcu- 
lated among them in order to obtain a median frame to use 
for a more complete star finding. The star list obtained in 
this way was then used as input for ALLFRAME, which 
performed simultaneous PSF fitting photometry on the 3 
V and 4 I frames. The whole procedure was repeated for 
each of the 8 chips. 

Calibration to the Johnson photometric system was 
performed by means of a set of Landolt (1992) standard 
fields observed the same night through the 8 WFI chips. 
The number of standard stars present in each of these 
fields has been recently increased by Stetson (2000), al- 
lowing us to measure ~ 40 standard stars per chip. A 
zero point and a color term were determined separately 
for each chip. The color terms were then averaged, and 
new zero points were calculated imposing the former to 
be fixed. Variations <0.1 mag have been found among the 
photometric zero points of different chips. 

Completeness estimates were determined for the WFI 
data in the same way described above for the SOFI ones. 
WFI data are more than 50% complete above I ~ 22.5, 
and V ~ 23. 



3. The Color Magnitude Diagram 

3.1. Near- 1 R CMD 

Figure [| shows the ncar-IR CMD of all the stars mea- 
sured in the bulge field. The jump in the number of stars 
at J = 16 is due to stars brighter than J = 16 having 
been sampled from the 8f3 x 8.'3 SOFI-LARGE field, while 
fainter stars have been measured with the larger angular 
resolution and deeper exposures of the SOFI-SMALL field, 
which mapped a 3'9 x 3'.9 area. 

The bulge HB red clump is visible at J ~ 14 and 
J — H = 0.6, partially merged into the RGB. The large 
magnitude spread of the HB is due to a combination of 
differential reddening, metallicity dispersion and depth ef- 
fect. These factors also cause the HB clump to merge ver- 
tically with the RGB bump, expected to be located ~ 0.7 
magnitudes fainter in J. The RGB bump (Iben 1968; Rood 
1972; Salaris, Cassisi & Weiss 2002) is predicted to be very 
populated in a high metallicity system as the bulge, and 
would have itself a large magnitude spread due to the same 
depth and reddening effects mentioned for the HB clump. 

The bulge turnoff is clearly visible at J ~ 17.5 and 
J — H ~ 0.4, while the almost vertical sequence depart- 
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Fig. 4. J, if (left) and J, A' s (right) CMD of the bulge 
stars imaged with SOFI. In both panels, stars brighter 
than J — 16 come from the 8.'3 x 8.'3 area mapped with 
short exposures and the SOFI-LARGE field, while fainter 
stars were measured in a 3.9 x 3.9 arcmin area observed 
with the SOFI-SMALL field. 



ing from near the bulge turnoff and extending upwards 
and bluewards is due to foreground main sequence stars 
belonging to the disk, widely dispersed along the line of 
sight (Ng & Bertelli 1996). 

Fig. ||a shows the CMD of the bulge field with super- 
imposed a 1 Gyr isochrone, for a solar metallicity popu- 
lation, adopting a distance modulus of (m — M)q = 14.47 
and a reddening of E(B — V) = 0.45 (see below). Clearly, 
a population of young stars physically located inside the 
bulge, would be definitely bluer than the vertical sequence 
of disk stars. The disk field, 30° away from the Galactic 
center, was used to statistically decontaminate from the 
foreground disk stars the bulge CMD shown in Fig. [|a. 
Small differences in the disk stellar population may be 
expected between the two lines of sight. Yet, the results 
have shown that the procedure is indeed quite effective. 
The CMD of the disk control field for an area of 24 square 
arcmin is shown in Fig. |E|b. The difference in reddening 
between the two fields was compensated by shifting the 
disk CMD shown in Fig. |b by 0.17 mag in J — H so as to 
match the color location of the blue disk sequences, and 
in magnitude by the corresponding Aj = 3.06E(J — H) 
extinction (Cardelli, Clayton & Mathis 1989). A region of 
the CMD free of bulge stars was then selected in order to 
normalize the number of disk stars observed in the disk 
field to the number of disk stars contaminating the bulge 
field. This region is indicated by the box in the upper left 
region of the CMD in Fig ||a. This normalization is correct 
only if no bulge star is present in the box. 

There are 91 stars in the box of the disk control field, 
which is 4.5 times less than in the corresponding box of 
the bulge SOFI-LARGE field. Therefore, 4.5 stars have 
been subtracted from the brighter region (J < 16) of the 



bulge CMD for each star seen in the disk CMD. Since the 
fainter part of the CMD (J > 16) was derived from the 
SOFI-SMALL field, the 4.5 scaling factor was divided by 
the ratio of the SOFI-LARGE to SMALL field area (4.6). 
Hence 0.98 stars have been subtracted from the fainter 
part of the CMD for each star in the disk control field. 

For each disk star in the disk CMD (Fig. ||b) we picked 
up the closest star in the bulge CMD (see below), and 
subtracted it according to the normalization factors given 
above, and to the slightly different completeness of the 
two fields. The distance on the CMD from a disk star to 
each bulge star was defined as: 



I d = y/[7 x A( J - H)f + AJ 2 . 

The color difference has been enhanced by a factor of 7 
because the color is much less sensitive than the magni- 
tude to physical differences (in the distance, reddeninig or 
mass) between a given disk star in the control field, and 
another disk star along the bulge line of sight. The result- 
ing, clean CMD of the bulge is shown in Fig. ||c, while 
the CMD of the stars statistically removed from the bulge 
CMD is shown in Fig. [|d. Also shown in this panel are 
the fiducial ridge line, helping the eye to identify the mean 
branches of the bulge CMD, and the brighter and fainter 
limits of our turnoff magnitude estimate (J = 17.65±0.2). 

The clumpy appearance of the RGB of Fig. ||d is due 
to the poor statistics in this region of the disk CMD since 
for each disk star we had to subtract 4.5 times more stars 
in the bulge field. Hence for each disk star a small "clus- 
ter" of close stars was subtracted from the bulge CMD. In 
the attempt to minimize this effect we actually subtracted 
every other three closest stars in the region J < 16. We 
expect a negligible effect on the 0.25 magnitude bins of 
the luminosity function discussed in the next Section. 

The 2MASS catalog was used to improve the statistics 
for the bright part of the CMD and of the luminosity func- 
tion. Note however, that the ~ 4" resolution of 2MASS 
is significantly worse than either the SOFI or WFI data. 
This complicates somewhat the cross identifications with 
the SOFI and WFI databases, and crowding effects are 
much more severe since two objects closer than ~ 4" are 
counted as one in the catalog. 

From the 2MASS point source catalog we extracted 
near-IR data for two large areas, 31 x 31 arcmin each, 
containing our bulge and disk field, respectively. A proce- 
dure identical to that described above for the SOFI data 
has been applied to the bulge CMD of Fig. ||a in order 
to statistically decontaminate bulge from disk stars, using 
the disk CMD shown in Fig. ||b. Note that since the disk 
and bulge field areas are now identical, the scaling factor 
for the decontamination (determined within a box similar 
to that shown in Fig. ||) is now very close to 1. 

The result is again very satisfactory (Fig. ||c,d) for 
stars brighter than J ~ 13. For fainter stars the decontam- 
ination is less effective because the magnitude and color 
spread of the disk and bulge CMD in Fig. ^|are quite differ- 
ent, due to the different drivers of the magnitude limits. 
Indeed, due to the low spatial resolution of the 2MASS 
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Fig. 5. a) The CMD of the bulge field. The solid line is the 1 Gyr isochrone for a solar metallicity population, b) 
CMD of the disk control field in the direction (I, b) — (30, 0). The region inside the box has been used to normalize the 
number of disk stars seen through the bulge line of sight, c) CMD of the bulge field as statistically decontaminated 
from the disk population. The horizontal lines locate the main sequence turnoff at J = 17.65 ± 0.2. The ridge line of 
the CMD is also shown, d) Stars that were subtracted from the bulge CMD in order to obtain the decontaminated 
CMD. 



data, the photometry is limited by the crowding in the 
bulge field, and by the background noise (including sky) 
in the disk. However, in the following only the brightest 
stars ( J<12) of the 2MASS decontaminated bulge CMD 
will be used. 



3.2. The Optical CMD 

Figure shows the (I, V - I) CMD of the 883,417 stars 
measured in the 34 x 33 arcmin field of WFI@2.2m. The 
eight panels in this figure correspond to the eight CCD 
chips of the WFI mosaic, mantaining their relative spatial 
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Fig. 6. CMDs of the bulge (panel a) and disk (b)) fields as constructed from the 2MASS photometry, c) The bulge 
CMD after statistical decontamination from disk stars; d) bona fide disk stars subtracted from a) in order to obtain 
c). 



position. The huge number of points per panel saturates 
the plot in the most populated areas, like the bulge main 
sequence (I > 18). The bulge turnoff is located around 
I = 18, while the main sequence of the foreground disk 
hits the bulge locus approximately midway between the 
turnoff and the base of the RGB. The HB red clump 
(merged with the RGB bump) is visible at / ~ 14.5 and 
V-I ~ 1.8, while the upper RGB (brighter than the HB) 
is extremely wide due to the bulge metallicity dispersion. 
Note, however, that the narrow sequence departing from 
the bulge HB and extending upward, almost vertically is 
due to the red clump of the disk stars, dispersed in mag- 
nitude as a result of their large spread in distance (and 
reddening). The bulge main sequence extends almost ver- 
tically in this plot, becoming very broad towards faint 



magnitudes for the combination of photometric error and 
plot saturation, but its increasing skewness towards red 
colors is in fact due to the presence of the faint extension 
of the disk main sequence. 

Although there is evidence for a bulge blue HB popula- 
tion in all the CMDs (with V-I ~ 0.5-=- 1 and I ~ 15H-18), 
note that chip#6 (panel f) contains the bulge globular 
cluster NGC 6558, already known to have a blue HB (Rich 
et al. 1998). The fact that the blue HB seen in this CMD 
belongs to the cluster becomes very clear when imposing a 
spatial selection around the cluster center. The other clus- 
ter features are not distinguishable from the bulge ones in 
this plot. A circular region of 2.4 arcmin radius, centered 
on NGC 6558, has been excluded from the following analy- 
sis. The radial trend of both cluster star counts and surface 
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Fig. 7. The optical CMD of the bulge from the WFI data. Each panel corresponds to each of the 8Kx8K CCD chips 
of the WFI camera. 



brightness flattens to the background level well inside this 
radius, ensuring a negligible cluster contamination outside 
2.4 arcmin. 

For a better visibility of the bright portion of the CMD, 
Fig. U displays only stars brighter than / = 16, this time 
combining data from all the 8 chips of WFI. The disk 
main sequence is very prominent on the left side of the 
diagram, and is paralleled some ~ 0.8 mag redder in V — I 
by the core helium burning clump sequence belonging to 
the same population. The bulge RGB and AGB occupy 
the right side of the diagram, and become very dispersed 
in their upper part due to the metallicity dispersion. The 
HB clump and the RGB bump are also indicated in the 
figure. 

4. The Metallicity Distribution of the Bulge 

From the combination of the optical (WFI) and near-IR 
(2MASS) decontaminated data, a (M K ,V-K) CMD was 
obtained including ~ 22, 000 stars with Mjc < 0.5, and the 
result is shown in Fig. ^|. A comparison with the fiducial 
line of the globular cluster NGC 6528 from Ferraro et al. 
(2000) immediately suggests that the mean metallicity of 
the bulge is virtually identical to that of this cluster, with 
a quite modest dispersion about this mean. This is demon- 




HB clump 
RGB bump 
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Fig. 8. The optical CMD for the bright stars combining 
the data in the whole WFI field. The location of the HB 
clump and the RGB bump is indicated. 
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Table 2. Globular Cluster used as RGB Templates 



Cluster 


[M/H] 


(m - M) 


E(B - V) 


NGC 6528 


-0.10 


14.37 


0.62 


NGC 6553 


-0.10 


13.46 


0.84 


NGC 104 


-0.57 


13.32 


0.05 


NGC 6171 


-0.88 


13.95 


0.31 


NGC 6121 


-1.06 


11.68 


0.35 


NGC 6809 


-1.59 


13.82 


0.10 


NGC 7099 


-1.84 


14.71 


0.01 


NGC 4590 


-1.90 


15.14 


0.04 


NGC 7078 


-1.92 


15.15 


0.09 




NGC 6538 " 



z I i I i i i I i i i I i i i I i i i I i i i I i i i I 

jfj 4 6 8 10 12 

(V-K) 

Fig. 9. Bulge CMD resulting from the combination of the 
optical WFI and near-IR 2MASS data. The line represents 
the fiducial locus of the globular cluster NGC 6528. The 
distance and reddening of the bulge field adopted for this 
comparison are (m - M) = 14.47 and E(B — V) = 0.45, 
respectively, while for NGC 6528 it is assumed (m— M)q = 
14.37 and E{B — V) = 0.62 (Ferraro et al. 2000). 



strated by the slope of the giant branch of the cluster and 
of that of the bulge field being the same. (Note that the 
metallicity affects the slope, while the reddening causes 
a solid shift of the RGB.) Therefore, the metallicity of 
NGC 6528 (and that of the twin cluster NGC 6553) has a 
special importance in connection with our attempt at de- 
termining the bulge metallicity distribution. In Fig. ^|, a 
distance modulus of (m-M) = 14.47±0.08 (7.8 kpc) was 
adopted for the bulge, according to the most recent deter- 
mination by McNamara et al. (2000) using RR-Lyrae and 
S Scuti variables from the OGLE survey. The reddening of 
our bulge field relative to NGC 6528 can then be estimated 
from Fig. ^, and we obtain E(B — V) = 0.45, needed to 



match the bulge and the cluster loci given the small dif- 
ference in the two distance moduli. For NGC 6528 the 
corresponding quantities are given in Table §. The large 
number of stars in the upper RGB of Fig. |, coupled with 
the high sensitivity of the (V — K) color to metallicity, 
and its very small sensitivity to a possible age dispersion, 
allows a determination of the bulge MD via the method 
described in Saviane et al. (2000). 

4.1. The Method 

The method is based on the construction of a family of hy- 
perbolas in the plane (Mx , V — K) suitable to represent 
a grid of upper RGBs from empirical template globular 
clusters, with known metallicities. Each hyperbola is rep- 



resented by the expression: 

M K = a + bx(V-K)+ {v _ C K) _ d (1) 

where the coefficients a, 6, c and d are quadratic functions 
of the metallicity: 

a = fci[M/H] 2 + fc 2 [M/H] + k 3 (2) 

fe = fc 4 [M/H] 2 + fc 5 [M/H]+fc 6 (3) 

c = fc 7 [M/H] 2 + fc 8 [M/H]+fc 9 (4) 

d = k w . (5) 



Finally, the inversion of equation (1) gives a value of the 
metallicity for each point in the (Mx, V — K) plane, hence 
for each bulge star. 

Although in principle spectroscopic determinations 
may be more accurate, the photometric approach has the 
advantage of relying on high S/N quantities (i.e., the 
magnitudes of the brightest stars) and on large statis- 
tics. However, it is entirely differential, i.e., it depends 
on the accuracy of the metallicity assigned to the tem- 
plates, which ultimately relies on spectroscopic determina- 
tions. The most complete and homogeneous grid of near- 
IR GC RGB templates presently available in the literature 
(Ferraro et al. 2000; see Table |]) includes 10 low redden- 
ing, nearby GCs with metallicity [M/H] between —1.92 
and —0.1. One cluster (NGC 6637) from the Ferraro et 
al. sample has been excluded from our calculations, for 
reasons explained below. 

With the exception of the two most metal rich clusters 
(NGC 6528 and NGC 6553), for the GCs in Table | we 
adopted the iron abundance [Fe/H] from Harris (1996), 
and corrected it to [M/H] adopting an a-element en- 
hancement [a/Fe]=0.3 for clusters with [Fe/H]< —1.0, 
and [a/Fe]=0.20 for more metal-rich clusters (Carney 
1996; Salaris & Cassisi 1996). As emphasized earlier, the 
case of NGC 6528 and NGC 6553 is particularly im- 
portant. A few stars in each of the two clusters have 
been observed at intermediate and high spectral res- 
olution, but different groups have obtained quite dis- 
crepant results. For NGC 6528, Carretta et al. (2001) 
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Fig. 10. The V, K RGBs of the template GCs published 
by Ferraro et al. (2000) (open squares), together with the 
analytical RGBs (dotted lines). The values of [M/H] for 
the analytical RGBs go from solar to [M/H] =—2.0 in steps 
of 0.2. 



and Coelho et al. (2001) report respectively [Fe/H]=+0.07 
and —0.5 (the latter value coming from low-resolution 
spectra). For [M/H] the same authors derive +0.17 and 
-0.25, respectively. For NGC 6553 Barbuy et al. (1999) 
give [Fe/H] = -0.55 and [M/H]=-0.08, while Cohen et al. 
(1999) report [Fc/H] = -0.16, and Origlia, Rich & Castro 
(2002) give [Fe/H]=-0.3, with [a/Fe]=+0.3. These dis- 
crepancies are uncomfortably large, and may hopefully 
disappear as soon as more high-resolution data are gath- 
ered at 8-10m class telescopes. For both clusters we finally 
adopt [M/H]=-0.1 (the value reported in Table 2) with 
the uncertainty of the resulting MD being dominated by 
the uncertainty of the metallicity of these two clusters. 

Figure [h] shows the resulting grid of RGB loci: open 
squares represent the fiducial points extracted from the 
empirical templates, while the RGBs from the analytical 
interpolation are shown as dotted lines. The coefficients 
of the Equations 4-7 are listed inside the figure. The an- 
alytical RGBs are shown for metallicities ranging from 
[M/H]=0 to [M/H]=-2.0, in steps of 0.2. All the empirical 
templates, with the exception of NGC 104 (third cluster 
from the right) follow very well the RGB shape trend of 
the analytical solution. We could find no obvious explana- 
tion for the apparent discrepancy of NGC 104: its CMD is 
quite populated in this region, and, being one of the best 
studied clusters, the adopted parameters (metallicity, red- 
dening and distance) are rather robust. The near-IR ob- 
servations of this cluster were done with the ESO/MPG 
2.2m telescope mounting IRAC-2, an early IR array de- 



Fig. 11. Upper panel: the CMD of the globular cluster 
NGC 6553 with overplotted the analytical RGB templates. 
The dispersion in distance typical of the bulge field at 
b = —6° and a reddening dispersion of AE(B — V) = 0.3 
have been added to the cluster CMD in order to simulate 
the observational biases affecting the derived bulge MD. 
Lower panel: MD obtained for NGC 6553. Shaded his- 
togram is the MD obtained if only the stars brighter than 
Mk = —4.5 are considered. This MD is well represented 
by a Gaussian centered in [M/H]=— 0.1 and a = 0.1. 



tector (Ferraro et al. 2000); perhaps the use of a more 
modern instrument may clarify the issue. 

As evident from Figure [l^, the sensitivity to metallic- 
ity is a strong function of the position in the CMD, hence 
of metallicity itself. The leverage increases with both lumi- 
nosity and color, being favored at high metallicities by the 
increasing TiO blanketing in the V band causing redder 
and redder V — K colors. 

In order to investigate the effect of distance and red- 
dening dispersion on the MD derived using the adopted 
method we simulated such effects on the (Mk,V — K) 
CMD of the bulge cluster NGC 6553 (Guarnieri et al. 
1998). To each star of NGC 6553 we associated a dis- 
tance modulus randomly extracted according to a bulge 
space density distribution of the form p = Const x r -3 7 
(Sellwood & Sanders 1988; Terndrup 1988), where r is the 
spatial distance from the Galactic center. The line of sight 
integration at the field position ((/, b) = (0,-6)) gives a 
1-tr dispersion of 0.13 mag in the distance modulus, the 
distribution being close to Gaussian in the core, but with 
somewhat broader wings. A top hat reddening distribu- 
tion of ±0.15 was also adopted. The latter is believed to 
be an upper limit to the observed bulge reddening varia- 
tion across our field, but has been chosen as a conserva- 
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tive assumption, in order to have an upper limit on the 
possible biases in the derived MD. Random extractions 
of both distance modulus and reddening were repeated 
6 times for each star in the original CMD by Guarnieri 
et al. (1998) in order to obtain a more populated CMD. 
Figure [ll] shows the resulting simulated CMD including 
distance and reddening dispersion effects (upper panel). 
The lower panel shows the MD obtained in the same way 
as that to be applied to the bulge, and reveals two inter- 
esting effects. First, if all the stars brighter than Mk = —3 
are considered, the resulting MD (dashed histogram) is ex- 
tremely broad with a long metal-poor tail. This is entirely 
due to the fact mentioned above: the separation between 
the RGBs corresponding to different metallicities becomes 
very small towards faint magnitudes, much smaller than 
the magnitude/color dispersion due to distance and red- 
dening dispersions. However, if the analysis is restricted 
to stars brighter than Mk = —4.5 then the derived MD 
shows a fairly sharp peak at [M/H]=— 0.1 (i.e., centered 
on the adopted metallicity for this cluster) , well fitted by a 
Gaussian distribution with a = 0.1 (FWHM=0.24). This 
exercise demonstrates that, as long as the analysis is con- 
fined to the stars brighter than Mk — —4.5: i) one should 
not expect systematic biases; ii) the spread introduced 
by the distance and reddening dispersions is of the order 
of ~ 0.1 dex in [M/H]; and Hi) contamination by AGB 
stars (which would introduce a bias artificially skewing 
the distribution towards the metal-poor side) is negligible 
at these magnitudes due to the exclusion of the more pop- 
ulated (and bluer) AGB "clump", located at Mk ~ —3. 
We therefore restrict to Mk < —4.5 the derivation of the 
bulge MD. 



4.2. The Resulting Metallicity Distribution 

Figure |lj shows the analytical RGBs overplotted to the 
bulge upper RGB. The dotted part of each RGB is the 
extrapolation above the theoretical location of the RGB 
tip (c.f., Section 7), hence only AGB stars are expected 
in this region. The metallicities of the analytical RGBs 
range from [M/H]=+0.2 to -1.8, in steps of 0.2. Note 
that the most metal-rich templates have metallicity of 
[M/H]=— 0.1, therefore the shape of the analytical RGBs 
for higher metallicity is the result of a small extrapolation. 
Stars bluer than V — K — 2.8 (vertical dashed line) were 
not considered, in order to exclude a few bluer objects that 
are most likely a disk residue due to imperfect statistical 
decontamination. Stars fainter than Mk = —4.5 (horizon- 
tal dashed line) were excluded as well, for the reasons ex- 
plained above. After these cuts, a sample of 503 stars has 
been used to derive the bulge MD. It is already clear from 
this figure that most of the bulge stars are located between 
the [M/H]~ 0.1 and [M/H]=-0.6 templates, with a peak 
at [M/H]=— 0.1 and very few stars more metal-poor than 
~ — 1. Formally, super-solar metallicity stars ([M/H]> 0) 
represent 29% of the sample. This should be regarded as 
an upper limit, due to the diffusion to higher [M/H] values 




(v-K)„ 

Fig. 12. The bulge V, K CMD compared with the ana- 
lytical RGB templates. The adopted bulge distance mod- 
ulus is (to — M)o — 14.47 while the mean reddening is 
E(B — V) = 0.45. Only the stars inside the dashed box 
where used for the MD determination. 



caused by the dispersion in distance, reddening, etc. (see 
Fig.0). 

Before proceeding to derive the bulge MD a few biases 
must be quantitatively evaluated and compensated for: 1) 
the rate at which stars leave the main sequence (called 
evolutionary flux in Renzini & Buzzoni 1986) has a slight 
metallicity dependence, 2) the rate of evolution along the 
RGB scales as Z-° M (Rood 1972), and 3) by adopting 
the Mk = —4.5 cut one samples (bolometrically) deeper 
on the RGB the higher the metallicities. 

The amplitude of these factors is shown in Fig. [l3| 
The upper panel shows how the evolutionary flux changes 
as a function on metallicity, relative to its value at Zq. 
The middle panel shows the combination of effect 2) and 
3) mentioned above: the time spent on the RGB, from 
Mk = —4.5 to the tip, normalized at solar value, slowly 
increases as a function of metallicity. Finally the lower 
panel shows the combination of the two factors above, 
which is the correction to apply to each metallicity bin of 
the MD to take into account the evolutionary effects. Since 
the evolutionary flux and the RGB lifetime have opposite 
behaviour with metallicity, the final correction factor is 
never larger than few percent, well within the intrinsic 
precision of our measurement. 

The resulting bulge MD for the 503 stars in Fig 
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is 



finally shown in the three panels of Fig. [l4| as a shaded his- 
togram. Also shown in Fig. [l4|, as thick "open" histograms, 
are the bulge MDs determined from various spectroscopic 
studies, normalized to the total number of stars in our 
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Fig. 13. Upper panel: the ratio of the evolutionary flux at 
any Z to the same quantity at solar Z . Middle panel: the 
time spent on the RGB from Mk = —4.5 to the tip, with 
respect to the value at Zq. Lower panel: the product of 
the two factors above, i.e., the correction factor that one 
should apply to the MD found from the CMD, in order to 
take into account these effects. 



sample. Note that all these spectroscopic determinations 
of the bulge MD refer to [Fe/H] abundance. Since the tem- 
plate RGBs are on the [M/H] scale, we subtracted from 
the bulge [M/H] distribution the same a-element enhance- 
ment (i.e., 0.2 for [Fe/H]> -1 and 0.3 for [Fe/H]< -1) 
that had been applied to the template GCs, in order to 
obtain a [Fe/H] distribution. Although the bulge high 
resolution studies do not permit to draw strong conclu- 
sions on the a-element enhancement, the similarity of the 
bulge CMD with that of bulge GCs would favor a similar 
chemical enrichment history. Certainly the new genera- 
tion multi-object spectrographs will allow us to improve 
our knowledge in this field. 



The comparisons shown in Fig. 14 shows that the 
present, "photometric" MD is broadly consistent with 
the spectroscopic ones by McWilliam & Rich (1994) and 
Ramirez et al. (2000), with just a somewhat less pro- 
nounced supersolar [Fe/H] tail in the photometric MD. 
(This characteristic is exacerbated in the comparison with 
the MD of Sadler et al. (1996), based on low resolu- 
tion data.) We caution, however, that the position of the 
high [Fe/H] cutoff is entirely dependent on the metallic- 
ity assigned to the two template clusters NGC 6528 and 
NGC 6553. For instance, the apparent discrepancy with, 
e.g., McWilliam & Rich (1994) at high metallicity would 
disappear if we would have assigned to NGC 6528 and 
NGC 6553 a metallicity ~ 0.2 dex higher. Indeed this 



Fig. 14. Comparison between the bulge MD derived from 
Fig. [12] (shaded histogram) and those derived with spec- 
troscopic surveys. 



would bring to coincidence with M&R the high metallicity 
cutoff, while stretching the distribution and reducing the 
excess at [Fe/H]~ —0.2. We believe that the homogeneous 
high-resolution studies of cluster and bulge field stars that 
will soon become available, e.g., with the forthcoming 
FLAMES multifibre spectrograph on VLT (Pasquini et 
al. 2000), will clarify this problem. 

On the other hand, we tend to consider as quite re- 
liable the sharpness of the cutoff at high [M/H]. Due to 
the TiO blanketing in the V band, the V — K color of 
RGB stars has a strong, rapidly increasing sensitivity to 
[M/H], and any high metallicity tail would have resulted 
in a dramatic broadening of the upper RGB, which in- 
stead is clearly absent in the CMD of the bulge shown in 
Fig. |9[ Moreover, the high [Fe/H] cutoff would be even 
sharper, had we underestimated the effect of the distance 
and reddening dispersions. 

4.3. Inferences for the Chemical Evolution of the Bulge 

The most straightforward approach is now to compare this 
empirical MD with the Simple (or one-zone, closed-box) 
Model of chemical evolution (Searle & Sargent 1972). Rich 
(1990) first made this comparison for the bulge and found 
the Simple Model to be a good fit to the abundance distri- 
bution from Rich (1988). The Simple Model assumes that 
no infall or outflow of metals has occurred during the star- 
forming phase of the system. The distribution follows the 
relation N(Z) — y~ 1 exp~ z ^ v , where the single parameter 
y, the yield, is given by y = (Z). Hence, the yield is the 
average metallicity of the stars in a one-zone system fol- 
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Fig. 15. Comparison between the metallicity distribution 
of bulge giants from the present study (histogram) and the 
prediction of the Simple Model {closed-box model) with 
the indicated value of the yield (y = (Z)). The model 
distribution has been scaled by the total area under the 
data histogram. There is a clear shortage of stars near zero 
metallicity, and an excess of stars near solar abundance, 
but the one-zone model appears to fit the abundance dis- 
tribution reasonably well. 



lowing the exhaustion of the gas (Hartwick 1976; Mould 
1984). If metals are lost from the system by outflow, the 
functional form of the distribution is unchanged, but the 
mean abundance (Z) is the yield (as would be calculated 
from supernova models) reduced by the appropriate out- 
flow parameter. 

Fig. [l5] shows a plot of the Simple Model distribution 
(for y — (Z) — 0.015) overlaying the data. Note that both 
panels show the same model and the same data, once us- 
ing [M/H] and once using Z/Zq as abscissa. As the Simple 
Model satisfies the requirements of a probability distribu- 
tion, the function is scaled by the area under the data 
histogram, so that the area covered by data and fit are 
identical. 

From the plots, the general shape of the abundance 
distribution is in fairly good agreement with the Simple 
Model. However, there are noticeable deviations: in the 
subsolar regime (Z < 0.3 Zq) less stars are observed than 
predicted, while with respect to the model there is an ex- 
cess of near-solar metallicity objects. This excess can be 
traced back to the fairly sharp clustering of the bulge 
stars around the RGB template at such metallicity (see 
again Fig. 7 and 10). This feature is most likely real, as 
many factors conjure to broaden the derived distribution 
(dispersion along the line of sight, differential reddening, 



Fig. 16. The bulge near-IR LF as obtained from the SOFI 
data. 

incomplete disk decontamination, etc.). For the same rea- 
sons, most likely real is the sharp cutoff at high metallic- 
ity, that may even be sharper given the mentioned effects 
that tend to smooth out any sharp feature in the MD. 
The moderate shortage of metal-poor stars compared to 
the Simple Model may flag the presence of a G- dwarf prob- 
lem, a point on which we return in Section 8. 

The sharp high-metallicity cutoff of the MD suggests 
that star formation did not proceed to complete gas con- 
sumption. If the bulge formed in a rapid and intense star- 
burst, one can imagine that much of the metals would 
have been produced in situ in core-collapse supernovae, 
eventually driving strong, metal-rich galactic winds. There 
are at least two independent arguments favoring this sce- 
nario, one direct and one indirect: a) high redshift galax- 
ies thought to be proto-bulges are observed to have strong 
metal enriched winds (Steidel et al. 1996), and b) in clus- 
ters of galaxies at least as much metal mass is out of galax- 
ies in the intracluster medium, as there is locked into stars 
inside galaxies, which is taken as evidence for early metal- 
enriched galactic winds (Renzini 1997, 2002). Therefore, 
a relatively sharp cutoff at the metal rich end, such as 
we observe, would result if the star formation was suffi- 
ciently violent as to eventually evacuate the gas from the 
proto-bulge before it was exhausted by star formation. 

5. The Luminosity Function of the Bulge 

The disk-decontaminated CMD shown in Fig. has been 
used to construct the bulge luminosity function (LF) from 
the RGB tip down to ~ 1 mag below the main sequence 
turnoff, in each of the three bands: J, H and K Sl with 
the results being shown in Fig. ufl. The number counts in 
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Fig. 17. The bulge LF from this work complemented with 
the one obtained in Paper II from NICMOS data. 

the SOFI-SMALL field (J > 16, J > 15.61) have been 
normalized to those in the SOFI-LARGE field, according 
to the ratio of their areas. 

For an easier comparison with the LFs of other ob- 
jects, or different bulge regions, the LFs presented here 
are always shown as a function of dereddened magnitude. 
The adopted average E(5 — V) is 0.45, as derived from the 
comparison with the fiducial loci of the CMD of NGC 6528 
(Fig. [|). The relations between the absorptions in differ- 
ent bands have been assumed to be Ay = 3.1 x E(B — V), 
Ai = 0.479 x A v , Aj = 0.282 x A v , A H = 0.190 x A v , 
and Ak = 0.114 X A v (Cardelli et al. 1989). 

The broad peak at J ~ 13.2 in Fig. |l6| is the HB 
clump, slightly "bimodal" because merged with the RGB 
bump, as already discussed in Section ^. The main se- 
quence turnoff is located at Jo = 17.05 ± 0.20, as deter- 
mined from Fig. |[ The steep decrease in the number of 
stars that could be expected just above the turnoff, is in 
fact smeared in this LF due to depth effect, differential 
reddening, metallicity dispersion, photometric errors and 
blending effects (see below). 

Figure [l7] shows the comparison, in J and H, be- 
tween the SOFI LF and the very deep NICMOS LF from 
Paper II. A simple normalization by the ratio of the field 
areas brings the NICMOS LF in perfect agreement with 
the SOFI one, based on a field more than 100 times larger. 
Note that the disk contribution to the NICMOS LF was 
subtracted using Kent (1992) model LF of disk and bulge 
(cf., Paper II for details). 

The upper panel of Fig. [l8| shows the comparison 
between the SOFI LF and the LF by Tiede, Frogel & 
Tcrndrup (1995; TFT) in the same bulge region, but on 



Fig. 18. Comparison with previously published bulge LFs. 
Counts from different sources have been scaled accord- 
ing to the ratios of the field area and to the different 
surface brightness, if referring to another bulge region. 
Histogram: this work; filled squares: data from Tiede, 
Frogel & Terndrup (1995); open circles, data from Frogel 
& Whitford (1987); filled triangles: data from the 2MASS 
sky survey. 

a smaller area (4056 arcsec 2 ). The disk contribution has 
been subtracted from the TFT star counts using the ratio 
between bulge and disk stars computed from Fig. ||a,b. 
The upper panel of Fig. [l8| also shows the LF obtained by 
Frogel & Whitford (1987; FW in the figure label) for the 
M giants in Baade's Window. The latter has been normal- 
ized both for the different area and the surface brightness 
difference between the two fields. 

The lower panel of Fig. [l8| shows the comparison with 
the counts from the 2MASS survey, also normalized only 
for the ratio of the field areas. The 2MASS counts plotted 
here were extracted from a region of 927 arcmin 2 , centered 
on the SOFI field. In the range where they are complete, 
i.e., for Jo < 15 the 2MASS counts agree perfectly well 
with the SOFI counts, and are consistent with TFT and 
FW counts. 

Finally, Fig. [l^ shows the comparison with the K-b&nd 
bulge LF from DePoy et al. (1993). The latter was ob- 
tained from the photometry of a 604 square arcmin fields 
towards Baade's Window, and included disk stars. The 
two LF where normalized according to the different field 
area and surface brightness, and the foreground disk con- 
tamination, estimated from the present data, has been 
subtracted from the Depoy LF. 

From the combination of the data from 2MASS, SOFI, 
and NICMOS, a composite LF for the Galactic bulge can 
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-0.184 


7.86 
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7.52 
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8.11 


0.144 
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0.418 


8.98 


0.027 


8.61 


0.504 


8.27 


0.469 


9.23 


0.621 


8.86 


0.526 


8.52 


0.418 


9.48 


0.469 


9.11 


0.536 


8.77 


0.481 


9.73 


0.547 


9.36 


0.594 


9.02 


0.557 


9.98 


0.594 


9.61 


0.805 


9.27 


0.776 


10.23 


0.566 


9.86 


0.794 


9.52 


0.800 


10.48 


0.811 


10.11 


0.962 


9.77 


0.904 


10.73 


1.007 


10.36 


1.062 


10.02 


0.926 


10.98 


1.027 


10.61 


1.138 


10.27 


1.144 


11.23 


1.136 


10.86 


1.149 


10.52 


1.095 


11.48 


1.128 


11.11 


1.210 


10.77 


1.260 


11.73 


1.279 


11.36 


1.352 


11.02 


1.250 


11.98 


1.352 


11.61 


1.383 


11.27 


1.428 


12.23 


1.415 


11.86 


1.447 


11.52 


1.416 


12.48 


1.633 


12.11 


1.553 


11.77 


1.500 


12.73 


1.623 


12.36 


1.733 


12.02 


1.568 


12.98 


1.857 


12.61 


2.003 


12.27 


1.784 


13.23 


2.068 


12.86 


2.075 


12.52 


2.006 


13.48 


2.013 


13.11 


2.090 


12.77 


2.085 


13.73 


2.143 


13.36 


2.114 


13.02 


2.121 


13.98 


2.114 


13.61 


2.100 


13.27 


2.188 


14.23 


1.959 


13.86 


1.959 


13.52 


1.991 


14.48 


2.033 


14.11 


1.982 


13.77 


1.996 


14.73 


1.944 


14.36 


1.892 


14.02 


1.978 


14.98 


2.068 


14.61 


2.072 


14.27 


1.919 


15.23 


2.140 


14.86 


2.111 


14.52 


2.033 


15.48 


2.188 


15.11 


2.083 


14.77 


2.090 


15.73 


2.412 


15.36 


2.487 


15.02 


2.179 


15.98 


2.756 


15.61 


2.665 


15.02 


2.187 


16.23 


2.790 


15.86 


2.818 


15.27 


2.422 


16.48 


2.957 


16.11 


2.918 


15.52 


2.740 


16.73 


3.130 


16.36 


3.071 


15.77 


2.813 


16.98 


3.228 


16.61 


3.191 


16.02 


2.916 


17.23 


3.376 


16.86 


3.320 


16.27 


3.080 


17.48 


3.518 


17.11 


3.478 


16.52 


3.218 


17.73 


3.666 


17.36 


3.598 


16.77 


3.285 


17.98 


3.612 


17.61 


3.676 


17.02 


3.504 


18.23 


3.591 


1 / .8o 


O TOT 

6. (of 


1 T 07 


aa 1 
O.OOl 




9 R7R 

O.Ui u 


18 11 


o.ouo 


1 7 ^9 


*3 7fi0 
O. i uu 


18 7^ 


O. i O^i 


1 8 


^ fton 


1 7 77 
li.fi 


*3 ft09 


1 8 Q8 


■3 81/1 
O.Ol4 


1 8 f!1 
lo.Dl 


^ Q77 

o.y it 


1 ft 09 


*3 ft9fi 


1 Q 


O.OD4 


1 8 8fi 
lo.oO 


/i nn 

4.Ulo 


1 ft 97 


*3 QO/1 


1 Q 48 


o.yui 


1 Q 1 1 




1 p:o 


■3 Q9fi 
o.yzu 


in 79 
iy. / o 


o.yoz 


1 Q ^fi 

iy .ou 




1 ft 77 
lo. i i 


■3 Q71 

o.y / 1 


1 Q Q8 




1 Q fil 

ly .01 


A OQ/1 


1 Q 09 


*3 QO*3 


20.23 


4.039 
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4.138 
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20.36 
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a Horizontal lines mark the boundaries between 2MASS, 
SOFI and NICMOS based data, from top to bottom. 
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Fig. 19. Comparison with the if-band LF by DePoy et al. 
(1993). Open symbols refer to counts affected by incom- 
pleteness at the faint end and by saturation at the bright 
end. 

be constructed, using the best data for each luminosity 
range. The result is reported in Table ||, that lists the star 
counts for J, H, and K s bands, while Table ?? lists the 
optical V and I LF. Note that NICMOS data are available 
only in the J and H bands. All the counts have been 
normalized to the area of Sf 3 x 8f 3 mapped by the SOFI- 
LARGE field: i.e., the 2MASS counts have been divided 
by 13.43, the SOFI-SMALL counts have been multiplied 
by 4.6, and the NICMOS counts have been multiplied by 
609. These scaling factors can be used to calculate the 
Poissonian errors associated with the counts in each bin. 
The numbers in Table [| have been corrected for both disk 
contamination and incompleteness, although the latter is 
only significant for the few faintest bins of the NICMOS 
LF, and it is always <50%. 

5.1. Simulated CMD and Theoretical LF 

In order to compare our observations to corresponding 
theoretical predictions we have developped a simulator 
which generates the CMD of a stellar population with 
a single age and a wide metallicity spread. In this way, 
we neglect the presence of an age spread, which is jus- 
tified since the location of the RGBs of (relatively) old 
stellar populations is much more sentitive to metallicity 
than to age. The code results from a development of the 
CMD simulator used in Greggio et al. (1998), adapted to 
describe single age Stellar Populations with a wide metal- 
licity spread. A thorough description of it can be found in 
Rejkuba (2002, PhD thesis). We give here a short report, 
and specify the ingredients used in our simulations. 

A Monte Carlo procedure is used to extract mass 
and metallicity of a simulated star, which gets then po- 
sitioned on the H-R diagram via interpolation among 
isochrones. The determined bolometric magnitude and ef- 
fective temperature are transformed into monochromatic 
magnitudes by interpolating within bolometric correction 
tables. Incompleteness and photometric errors, as mea- 
sured on the real frames are applied. In particular this 
includes the brightening effect due to blending that was 
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Table 4. Optical bulge luminosity function. ~ 6 
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discussed in Section 2.1. The procedure is iterated until 
the number of objects observed in some CMD region is 
reached. 

The isochrones data base consists of Cassisi and Salaris 
(1997) models, implemented with Bono et al. (1997) , plus 
some additional models explicitly computed for this ap- 
plication with the code described in Cassisi and Salaris 
(1997) and Cassisi, Salaris and Bono (2002). The metal- 
licity range goes from Z=0.0001 to Z=0.04, and the he- 
lium abundance varies in locksteps with Z following Y ~ 
0.23 + 2. r oZ. 



0.5 1 

(J-K)„ 



0.5 1 
(J-K)„ 



Fig. 20. Left panel: simulated CMD for a 13 Gyr old pop- 
ulation with the MD determined in Section 4. Right panel: 
the bulge CMD for the SOFI-LARGE field. 



Theoretical bolometric corrections (BCs) have been 
obtained by convolving the model atmospheres computed 
by Castelli et al. (1997) with the Landolt V and / and the 
SOFI J, H, K s filter passbands and fixing the zero point 
such that BCv = —0.07 for the model of the Sun, and all 
the colors are zero for the model of Vega. However, the 
brightest stars in our sample are cooler than T ~ 4000K, 
limit below which current model atmospheres are known 
to fail to reproduce the observed spectra, due to the inap- 
propriate treatment of molecules. Therefore, the empirical 
BCs by Montegriffo et al. (1998) where used instead of the 
theoretical ones for temperatures T<4000 K. The depen- 
dence of the BCs on the metallicity parameter is taken 
into account in the simulations. 

The results of the artificial stars experiments described 
in Section 2 have been used to assign the detection prob- 
ability and the photometric error (i.e. the difference be- 
tween input and output magnitude) to the simulated stars. 
The synthetic CMD has thus the same observational bi- 
ases as that obtained from the measured frames. 

A number of simulations have been computed for ages 
ranging from 8 to 15 Gyr, and by adopting: 

i) an IMF slope of —1.33, as found in Paper II 

ii) the metallicity distribution derived in Section 4 

in) an average distance modulus of (m — M) = 14.47 and 
reddening of E(B — V) = 0.42Q Each obtained synthetic 



CMD has been further dispersed for depth ((p(r) oc 



„-3.7 



and variable reddening (AE(B — V) = ±0.15) effects. 

Figure |o| (left panel) shows the synthetic CMD for a 
13 Gyr old population with the characteristics described 



4 This value is 0.03 magnitudes smaller than the one adopted 
in Section 4, which was derived from the comparison of the 
mean locus of the (M V ,V - K) bulge CMD with that of 
NGC 6528. The lower value of E(B - V) = 0.42 is instead 
required for a best fit with the adopted isochrones, and there- 
fore for a consistent comparison with the synthetic CMD. 
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Fig. 21. The complete bulge LF in the three near-IR 
bands (dots), as resulting from the combination of the 
2MASS, SOFI and NICMOS data (from Table g). A the- 
oretical LF for an age of 13 Gyr (solid line) and 8 Gyr 
(dotted line) is shown as a solid line for comparison. 

above, and a comparison with the bulge near-IR CMD 
observed with SOFI is shown in the right panel. 

The simulation code was run until the number of stars 
with Mj < 1.2 reached the number of stars sampled by 
SOFI-LARGE in this CMD region. Then it was run again 
until extracting as many stars with Mj > 1.2 as in the 
SOFI-SMALL field. 

The main purpose of the comparison in Fig. is to 
show how much the simulation code is able to reproduce 
the observed CMD, including the observational biases in- 
troduced by the dispersions in distance, reddening and 
metallicity, as well as blending. 

Figure ^ shows that indeed the main features of the 
observed CMD are well reproduced by the simulator. 
Noticeable exceptions are the morphology of the HB 
clump, and of the lower RGB. In fact, the observed HB is 
significantly less defined than the simulated one, and also 
the color width of the RGB is underestimated by about 
a factor of two in J — K . In principle, this mismatch can 
be caused by several effect, although none of them, alone, 
seems convincing to us. Larger spread in the observed 
CMD can be due to an underestimate of one of the 
following effects: 

i) Distance spread: the bulge density law could be flat- 
ter than p oc r~ 3 ' 7 , and/or the bulge being a bar may in- 
crease the distance dispersion along the line of sight. Note 
that adopting a larger distance dispersion would smear the 
HB clump, making is more similar to the observed one, 



but would not have any appreciable effect on the spread 
in color of the RGB. Hence, a larger distance spread alone 
would not be sufficient to make the simulated CMD iden- 
tical to the observed one. 

ii) Differential reddening. Larger values of this param- 
eter would broaden the RGB of the simulated CMD, but 
also would cause the HB to appear tilted along the red- 
dening line, which is not seen in the observed CMD of the 
bulge, while it is prominent in e.g., NGC 6553 (Zoccali et 
al. 2001a). Also, the differential reddening needed to ac- 
count for the spread in color seen in the near-IR CMD 
would imply a too large spread in the optical CMD. 
An attempt was made to correct for differential redden- 
ing following the method used for NGC 6553 (Zoccali et 
al. 2001a), but failed due to the small scale of the red- 
dening variations across the field. Indeed, in the case of 
NGC 6553, one finds reddening variations of the order of 
SE(V - I) ~ 0.1 on scales of only - 20". 

Hi) Photometric errors. Larger SOFI-LARGE photo- 
metric errors than adopted in the simulation would cer- 
tainly help smearing the simulated RGB and HB respec- 
tively in color and magnitude, but can hardly solve both 
problems at once. To smear enough the simulated HB the 
photometric error should be ~ 0.3 — 0.4 mag, which would 
produce a too broad RGB compared to the observed CMD 
(see Fig. H). 

iv) Problems in the theoretical models. The RGB tem- 
peratures (hence colors) are strongly dependent on the 
mixing-length parameter, which needs to be empirically 
calibrated, and there are no perfect calibrators. The mod- 
els used in the present simulations were calibrated by fit- 
ting the RGB of template globular clusters from Frogel et 
al. (1983) for [Fe/H]< - 0.7, and at solar metallicity by 
demanding to the solar model to have the solar radius. 
Ideally, it would have been preferable to have a homoge- 
neous calibration over the whole metallicity range, but the 
solar metallicity clusters would have introduced in the cal- 
ibration the uncertainly in their reddening. As a net result, 
the adopted mixing length is some 10% larger at [Fe/H]=0 
than at [Fe/H]<— 0.7, which has the effect of compressing 
somewhat the relative spacing of the RGBs as a func- 
tion of metallicity. This effect goes indeed in the direction 
of reducing the color dispersion in the simulated RGBs. 
Moreover, the simulation is also affected by the uncer- 
tainty in the bolometric corrections and color-temperature 
transformations. All in all, Salaris, Cassisi & Weiss (2002) 
estimate an uncertainty of ~ 0.10 — 0.15 mag in the range 
of optical colors spanned by theoretical RGBs, which can 
also be taken as indicative of the uncertainty in the range 
of the near-IR colors. 

While the origin of the dicrepancy of Fig. is still 
partly unclear, we proceed to the construction of the theo- 
retical LF from the simulated CMD, keeping in mind that 
the region around the HB and lower RGB is presently 
not well reproduced by our simulations. The same code 
with the same inputs was then used to generate a much 
larger number of stars compared to the simulation shown 
in Fig. E(X in order to construct a smooth LF from the up- 
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per RGB down to the limit of the NICMOS photometry 
and the result is shown in Fig. The simulation includes 
photometric error and blending effects, and is meant to 
match the observed LF after correction for incomplete- 
ness. 

Figure |l] finally shows the comparison between the 
observed J, H , K s LFs (dots) and the theoretical LFs gen- 
erated from the simulated CMD described above (lines). 
The observed LFs result from the combination of all the 
available data, namely 2MASS + SOFI + NICMOS all 
scaled to the SOFI-LARGE area, as in Tables ||. The simu- 
lated LFs refer to 8 and 13 Gyr old populations. Observed 
and theoretical LFs were normalized to the total number 
of stars with MAGo < 12 in all bands. There is over- 
all agreement between the theoretical and the observed 
LFs, though the HB clump + RGB bump appears much 
sharper in the simulation than in the observed LF, as ex- 
pected from Fig. 

This comparison also shows that formally the 8 Gyr 
isochrone gives a slightly better fit to the data above 
the turnoff, which is the part of the LF sensitive to age. 
However, the same effect causing the smearing of the HB 
would have also made shallower the drop off of the lumi- 
nosity function just above turnoff, hence making the bulge 
population to appear younger than it is when comparing 
to the simulated LFs. Hence, we believe that before hav- 
ing identified the origin of the additional RGB and HB 
dispersion the comparison of simulated and observed LFs 
cannot set more stringent constraints on the age of the 
bulge stars other than being ~ 10 Gyr. 

5.2. The Bulge Spectral Energy Distribution 

Integrating the LFs in Tables || and [|, one can determine 
the total luminosity of the sampled bulge stellar popu- 
lation in each band. This has been performed using the 
whole database (NICMOS, SOFI, 2MASS for the infrared 
and WFI for the optical), then normalizing the results 
to the area of the SOFI-LARGE field. For the J and H 
bands the NICMOS data allow to include the contribution 
of all stars down to ~ O.15M . The A"-band contribution 
of stars fainter than the SOFI limit has been estimated 
from the theoretical M/L ratio of lower main sequence 
stars and adopting —1.33 for the slope of the IMF, and 
found to be just a fraction ~ 1CP 5 of the total AT-band lu- 
minosity. This result is not surprising given the flat bulge 
IMF. For the same reason, we integrated the V and / LFs 
of Table H safely neglecting the contribution of the lower 
MS stars. The resulting sampled luminosities are: 

L v = 1.00 x 10 5 i v ,o 
U = 1-13 x 10 5 £i, o 
Lj = 2.84 x 10 5 L JiQ 
L H = 3.67 x 10 5 L H ,o 
L K = 4.20 x 10 5 L K ,o, 

where the solar magnitudes have been taken from Allen 
(2000). To derive monochromatic fluxes, the integrated 
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Fig. 22. The spectral energy distribution of the bulge 
summing the contribution of all individual star (filled cir- 
cles). The error bars here show what would be the sys- 
tematic displacement of all the points if an uncertainty 
of E{B — V) = ±0.1 is allowed. Superimposed are the 
spectral energy distributions of template galaxies from 
Mannucci et al. (2001). 

magnitudes in each band have been transformed to AB 
magnitudes, and the monochromatic flux has been cal- 
cualted as logF,, = -0.4 m AB + 19.44 (Oke & Gunn 1983), 
with F\ — F p c/\ 2 . Figure ^2| shows the resulting spectral 
energy distribution (SED) of the galactic bulge, compared 
to the template spectra of elliptical, Sa and Sc galaxies 
from Mannucci et al. (2001). The Galactic bulge follows 
quite well the Sc template, which is not surprising given 
that the Milky Way is indeed an Sbc galaxy. For the op- 
tical part of the spectrum, the Mannucci et al. templates 
rely on those by Kinney et al. (1996), which in the case of 
Sc galaxies come from the integrated spectra of a region 
10" x 20" wide at the center of the two galaxies NGC 2403 
and NGC 598, whose radii are 10.6 and 35.2 arcmin re- 
spectively, i.e., the sampled regions are located well inside 
the bulge of these galaxies. Visual inspection of a DSS im- 
age of NGC 2403 suggests that the aperture used to derive 
the spectrum of this galaxy almost certainly includes a few 
site of active star formation. Their effect is however mini- 
mal in the V and / bands, while dominating the spectrum 
only in the IUE ultraviolet. 

On the other hand, the average value of the Mg2 in- 
dex of the ellipticals in the Kinney et al. sample is 0.314, 
which compares to Mg2=0.23 (Puzia et al. 2002) for the 
Galactic bulge. Hence, our bulge is quite metal poor com- 
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pared to giant ellipticals, which accounts for its bluer SED. 
Actually, the computation of the theoretical Mg2 index 
for a composite stellar population with the Bulge MD as 
determined here yields values as low as 0.16 — 0.19 for 
ages ranging from 10 to 15 Gyr (Maraston et al. 2002). 
These values, appropriate for simple stellar population 
with [Fe/H]~ —0.5, result from the effect of the well popu- 
lated subsolar metallicity component on the MD (Greggio 
1997). At optical wavelenghts, this component provides 
more relative flux than the high Z one, which results in 
both a lower composite Mg2 index, and in a bluer SED. 

5.3. The Sampled Luminosity-Star Number Connection 

The number of stars with mass in the range 0.15 to 1 M Q 
in the observed SOFI field can be obtained by integrating 
the bulge IMF with the appropriate scale factor A: 

Alfred = A [ M-^dM, (6) 

JO. 15 

with A tx Lt (Renzini 1998), where Lt is the bolometric 
luminosity sampled by the SOFI field. For given age, the 
ratio A/Lt is a weak function of both metallicity and 
IMF. Its dependence on age is relatively stronger, in fact, 
from the models by Maraston (2002), at [M/H]=-0.2 and 
for IMF slope of -1.33, A/L T ranges from 0.82 to 1.1 for 
an age varying from 10 to 13 Gyr. 

The total bolometric luminosity sampled by the SOFI- 
LARGE field can be obtained applying the appropriate 
bolometric corrections to individual stars in the sample. 
To evaluate that, we have run the simulation code to ob- 
tain the bolometric and monocromatic luminosities of a 
stellar population with the observed metallicity distribu- 
tion, obtaining: Lt = 1.39Ly = 1.16L/ = 0.70Lj = 
0.53Lh = 0.51Lk with the coefficients changing by at 
most 3% when varying the age from 10 to 13 Gyr. No 
observational errors were applied in this run. 

Using the five monochromatic luminosities with the 
above relation and averaging the results we obtain a total 
luminosity of L T = 177,000L o , hence Apred = 380,000 
and 510, 000, respectively for 10 and 13 Gyr. These num- 
bers compare with the 424,000±81, 000 stars in the SOFI- 
LARGE field, the error in the latter being dominated by 
the Poisson noise in the number of objects in the NICMOS 
field. Predicted and observed numbers are quite consistent 
(see also Paper I) , since these theoretical estimates are ex- 
pected to be accurate to within ~ 10%. 

5.4. The M/L Ratios 

Integrating the IMF from 0.15 to 100 Mq, and adopting 
the prescription in Paper II for the mass of the stellar 
remnants (those above ~ 1M Q ) one derives the total stel- 
lar mass in the bulge SOFI-LARGE sample. Using the 
luminosities in the various bands given above, one then 
determines the corresponding M/L ratios. This gives: 

M/L v = 3.67 
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Fig. 23. The bulge and the NGC 6528 CMDs are shown 
side by side. The magnitude difference between the HB 
and the turnoff of in the bulge CMD (left) is compared 
with the same quantity for NGC 6528 (right). The CMD 
for NGC 6528, originally in the NICMOS instrumental 
(mno,mno-mi6o) plane, has been shifted both in magni- 
tude and in color in order to match the location of the 
bulge HB. 

M/Li = 3.25 
M/L j = 1.28 
M/L H = 1-00 
M/L K = 0.87, 

with M/Lk in very nice agreement with the dynamical 
value, M/L K ~ 1 (Kent 1992). 

6. The Age of the Bulge 

Taking advantage of the sharp turnoff region of the decon- 
taminated near-IR CMD derived in Section 3 we proceed 
to estimate the age of the bulge stellar population. As in 
paper I, we adopt a differential procedure, comparing the 
luminosity difference between the HB clump and the MS 
turnoff of the bulge to that of a globular cluster of similar 
metallicity. 

Fig. ^3] shows the comparison between the bulge CMD 
and that of the cluster NGC 6528, whose metallicity is 
close to the average of the bulge. The near-IR CMD of 
NGC 6528 is based on the NICMOS photometry obtained 
by Ortolani et al. (2001). The magnitude difference be- 
tween the HB clump and the turnoff is virtually identical 
in the two diagrams, as emphasized by the two horizontal 
lines. 

This ensures that the difference between the age of 
the cluster and the mean age of the bulge cannot exceed 
~ 20% (thanks to the rule of thumb according to which 
5age/age~ S(AM$g) (Renzini 1991). 

This confirms and reinforces the conclusion in Paper I 
that the bulk of the bulge population and the clusters 
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Fig. 24. 10 Gyr isochrones (Cassisi & Salaris 1997) for 
the two extremes of the bulge MD are overplotted on the 
CMD. 

NGC 6528 and NGC 6553 are coeval. The absence of 
any appreciable extension of the bulge main sequence be- 
yond the obvious turnoff makes it clear that no trace of 
an intermediate-age population is detectable in the bulge 
CMD. 

The proper motion decontaminated and differential 
reddening corrected CMD of NGC 6553 (Zoccali et al. 
2001a) and NGC 6528 (Feltzing et al. 2002) confirm the 
results of Paper I that the HB to TO magnitude differ- 
ence in these two bulge clusters is virtually identical to 
that of the inner halo clusters NGC 104 (47 Tuc). Fig. |3| 
now shows that this magnitude difference is essentially 
identical also for the bulge field population, strengthen- 
ing the case for the bulk of the whole population of the 
Galactic spheroid (i.e. bulge and halo) being essentially 
coeval, though an age difference of ~ 20% (~ 2 — 3 Gyr) 
either way cannot be excluded. 

One aspect of the cluster to bulge CMD comparison 
still deserves some attention. Indeed, the bulge popula- 
tion is affected by dispersion in both distance and metal- 
licity, while the cluster stars are chemically homogeneous 
and all at the same distance (although affected by some 
differential reddening). In Fig. [24| two 10 Gyr isochrones 
spanning the full metallicity range of the bulge are over- 
plotted to the bulge CMD, assuming the same distance 
and reddening for both of them. This illustrates that the 
wider dispersion affecting the bulge CMD (compared to 
the HST/NICMOS CMD of NGC 6528) can be well ac- 
counted by the bulge metallicity dispersion, also taking 
into account the ~ 0.13 mag 1-cr dispersion due to the 
distance distribution along the line of sight. 




J-H J-H 

Fig. 25. Comparison of the bulge CMD with two younger 
isochrones of 3 (left panel) and 5 Gyr (right panel). Two 
models are plotted in each panel, both referring to the 
same age. The reddest curve in each panel is for so- 
lar metallicity, while the one on the blue side is for 
[M/H]=-1.3. 



According to recent attempts at determining the rela- 
tive ages of Galactic globular clusters the bulk of clusters 
are coeval within a ±1.5 Gyr uncertainty, with only the 
most metal rich ones in the sample appearing to be slightly 
younger than the others (Rosenberg et al. 1999; Salaris 
& Weiss 2002). However, these studies do not extend to 
the high-metallicity clusters of the bulge. For example, 
Rosenberg et al. assign to 47 Tuc an age 1.2 ± 1.2 Gyr 
"younger" than that the bulk of the halo globular clus- 
ters. Salaris & Weiss (2002) assign to the same cluster an 
age of 10.7 ± 1.0 Gyr, compared to 11.7 ± 0.8 Gyr for the 
prototypical metal poor cluster NGC 7078 (M15). On the 
other hand, Ortolani et al. (2001) date NGC 6528 at 13±3 
Gyr from the value of AJ™. 

It is clear that, within the uncertainties of the currently 
available data and dating methods, no appreciable age 
difference has been unambiguously detected between the 
bulk of bulge field stars and the globular clusters of either 
the bulge or the halo. On the other hand, the absolute 
age of the clusters remain more uncertain than the formal 
error bars sometime quoted by individual authors. Just 
to mention one example, the age of the globular cluster 
47 Tuc has been recently estimated to be 12.5±2 (Carretta 
et al. 2000), 13±2.5 Gyr (Zoccali et al. 2001b), and 10.7± 
1.0 Gyr (Salaris & Weiss 2002), the difference being partly 
due to a difference in the cluster distance and partly to 
the use of different sets of models. 

Significantly younger ages can be excluded, as shown 
in Fig. |25|, where 3 and 5 Gyr isochrones of both solar 
and [M/H]=— 1.3 metallicity are overplotted on the bulge 
CMD. After the submission of this paper, during the ref- 
ereing process, we became aware of the paper by Cole & 
Weinberg (2002), in which the Authors argue that the bulk 
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Fig. 26. The disk decontaminated optical CMD, together 
with the 3 Gyr isochrone for solar metallicity. The box 
shows the region where the stars of a hypotetical 3 Gyr 
old "bar" population, spread along the line of sight, would 
be expected. 



of the stellar population of the Galactic "bar" formed less 
than 6 Gyr ago, with an age of ~ 3 Gyr being favored. 
As they state "the main sequence turnoff of a 3 Gyr old 
population should be readily traceable along the Galactic 
bar from V ~ 17 at the near end to V « 19 at the far 
end" . Note that the Galactic component called "bar" in 
Cole & Weinberg (2002) has a mass of 2 x 10 10 M Q and 
therefore is not a minor component, but rather the whole 
population of the system called here "the bulge" . As evi- 
dent from Fig. ^6|, no such intermediate age population is 
actually detected in the present data. 

The region in the CMD just above the main sequence 
turnoff is so devoid of stars that very few, if any, blue 
stragglers stars (BSS) may be present in the field (see, 
e.g., Fig.^4|). Among Galactic globular clusters, Ferraro et 
al. (1995) estimate an average frequency of ~ 1 BSS every 
1O 3 L of bolometric light of the parent cluster, but with 
very large cluster to cluster variations that are not merely 
statistical fluctuations. Scaling from the SOFI-LARGE 
field, the SOFI-SMALL field samples - 177,000/4.6 ~ 
38, 000 L Q , and one would expect to recover ~ 38 BSSs, if 
the bulge has the same BSS frequency as the average glob- 
ular cluster. Clearly it has not. The bulge is far less pro- 
ductive of BSSs than a typical globular cluster, indicative 
that the cluster environment favors the formation of bi- 
naries with the right separation for producing BSSs. Most 
likely this is due to the dynamical processes that are ger- 
maine to the clusters. 



Fig. 27. The bulge RGB tip in different bands. The solid 
line is the theoretical tip location according to Salaris & 
Cassisi (1998), while the dashed line is the location of the 
tip empirically found by Ferraro et al. (2000). Big dots 
are the (20) stars that lie above the theoretical RGB tip 
in all the bands, while open circles are stars that lie above 
the tip in K but not in /. Crosses are blends even in the 
optical images. 

7. The RGB Tip and Above 

The tip of the RGB is of special interest for three main 
reasons: i) any AGB star brighter than the RGB tip is 
a candidate intermediate age star, hence may signal the 
presence of an intermediate age population (e.g. Iben & 
Renzini 1983), ii) the RGB tip itself is directly used as a 
distance indicator (Lee, Freedman & Madore 1993; Salaris 
& Cassisi 1997; Sakai, Madore, & Freedman 1999), and iii) 
in any galaxy the luminosity function near the tip of the 
RGB has a direct impact on the distance determination 
using the surface brightness fluctuation method (Tonry & 
Schneider 1988). 

Theoretical models indicate that the bolometric lumi- 
nosity of the RGB tip (corresponding to the helium flash 
in the core) increases with metallicity by ~ 0.25 mag 
per dex (Rood 1972; Sweigart & Gross 1978). However, 
due to increasing blanketing, in the / band the RGB tip 
magnitude appears to be independent of metallicity for 
[Fe/H]< — 0.7, i.e., in the range covered by globular clus- 
ters in the Galactic halo (Lee et al. 1993). The present 
data on the Galactic bulge offer an opportunity to ex- 
plore the behavior of the RGB tip luminosity at higher 
metallicities, up to solar and beyond. 

Fig. p7| a displays the (K, V — K) diagram of the 
combined WFI+2MASS field, for the upper RGB+AGB 
(K < 10). Overplotted, as a solid line, is the location 
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of the theoretical RGB tip for metallicities in the range 
from [M/H] = -2.35 to [M/H]=-0.28 (extrapolated here to 
higher metallicity) . For each metallicity, the RGB tip bolo- 
metric luminosities from the theoretical models of Salaris 
& Cassisi (1998) have been converted to K magnitudes 
and V — K colors using the empirical bolometric correc- 
tions from Montegriffo et al. (1998) and the analytical 
template RGBs described in Section 4. The empirical lo- 
cation of the RGB tip in the template globular clusters 
is also shown as a dotted line (from Ferraro et al. 2000), 
with the small offset between the theoretical and empir- 
ical RGB tips being probably due to poor statistics near 
the RGB tip in the template globular clusters. 

In Fig. [27], 25 stars lie above the nominal RGB tip in 
K, and therefore are candidate AGB stars. Three more 
stars with very red colors (J — K > 2.5) lie outside the 
limits of these plots. Two of them are included in the IRAS 
catalogue (IRAS18061-3140 and IRAS18069-3153) as OH- 
IR stars. The third one, missing in our optical CMD due 
to incompleteness (V mag for the other two stars is 20.7 
and 23.7, respectively) is likely to be of the same nature. 
Two stars in the same CMD region are also present in the 
disk CMD from 2MASS, therefore suggesting that they are 
not tracers of a younger bulge population. As well known, 
in metal poor globular clusters ([Fc/H]< — 1) no AGB 
stars brighter than RGB tip have been found. However, 
a few stars brighter than the RGB tip do exists in more 
metal rich globulars, with many of them being long pe- 
riod variables (LPV, Frogel & Elias 1988; Guarnieri et al. 
1997). The frequency of LPVs in metal rich globulars has 
been evaluated by Frogel & Whitelock (1998). From their 
Table 1, one derives an LPV frequency of ~ 4/(10 6 L^). 
The luminosity sampled by the combined WFI-2MASS 
field is 13.43 times larger than that sampled by the SOFI- 
LARGE field, i.e. ~ 5.7 x 10 6 L|, hence one would expect 
23 LPVs, while 28 stars are found above the RGB tip. 
However, an inspection of their optical counterparts in the 
WFI frames reveals that 6 of them, shown as crosses in 
Fig. are blended with other dimmer stars, and their in- 
trinsic iVband luminosity will be somewhat dimmer than 
estimated on 2MASS data. Another 5 stars are definitely 
below the RGB tip in the optical CMD (Fig. ||b ) , and are 
plotted as open circles in the four panels of Fig. E71 Their 
location may result either from the resolution of 2MASS 
being much worse (~ 4 arcsec) than that of WFI (<1 arc- 
sec) and therefore these stars may have unresolved com- 
panions in near-IR, that make them artificially brighter. 
Alternatively they might be LPVs having been caught 
near minimum light at the time of WFI observations. 
Finally, our simulations show that ~ 5 RGB stars are 
expected to be found above the nominal RGB tip due to 
the depth effect plus differential reddening. On the other 
hand, other LPVs may have been caught below the RGB, 
near minimum light. All in all we conclude that the num- 
ber of stars brighter than the RGB tip is within the expec- 
tations as derived from the frequency of LPVs in globular 
clusters spanning nearly the same metallicity range of the 
bulge. Of course, not all stars brighter than the RGB tip 
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Fig. 28. The absolute bolometric luminosity of upper 
RGB and AGB stars. The theoretical RGB tip luminos- 
ity is also shown as a solid line (the same as in Fig. p7|) . 
The bulge AGB stars appear to extend ~ 1 bolometric 
magnitude above the RGB tip. 



in Fig. ^ may be LPVs, hence the present findings do not 
contradict the result of Frogel & Whitelock (1998), who 
argue for a deficit of LPVs in the bulge field compared to 
metal rich globulars. In any event, the bright AGB stars 
seen in the WFI-2MASS field do not favor the presence 
of a stellar population in the galactic bulge appreciably 
younger than the old, metal rich globular clusters, thus 
reinforcing the conclusion we have drawn from from the 
turnoff region. 

The absolute bolometric magnitudes for the stars on 
the upper RGB and AGB have also been derived using the 
empirical bolometric corrections from Montegriffo et al. 
(1998) and the distance and reddening adopted in Section 
4. The result is displayed in Fig. [28], showing that the bulge 
AGB stars extend ~ 1 bolometric magnitude above the 
RGB tip, reaching Afboi — —5.0. This is just marginally 
brighter (by ~ 0.2 magnitudes) than the brightest known 
LPV member of a metal rich globular cluster, i.e. V4 in 
NGC 6553 (Guarnieri et al. 1997). 

Concerning the use of the RGB tip as a standard can- 
dle, Fig. |27j b shows that the I band magnitude of the tip 
drops by almost two magnitudes, from the lowest to the 
highest metallicity spanned by bulge stars. Clearly the I 
band luminosity of the tip cannot be used as a standard 
candle above [Fe/H]~ —0.7. Also in the J band the tip 



luminosity appears to decrease slightly (Fig. 27:), while 
it increases slightly in the K band (Fig. |27]a). In the H 
band instead the tip luminosity appears to be fairly con- 
stant, suggesting the use of this band for deriving RGB 
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tip distances of stellar populations in the metallicity range 
— 1.0<[Fe/H]<0.0. For even higher metallicities the K 
band may become more appropriate, as indicated by the 
theoretical RGB tip luminosity shown in Fig. |27|a. 

Fig. ^7] suggests also quite obvious considerations con- 
cerning the surface brightness fluctuation method. Clearly, 
the fluctuation magnitude will depend both on the partic- 
ular band used, and on the specific metallicity distribution 
in each galaxy. 

8. Summary and Conclusions 

Optical and near-IR CMDs for the Galactic bulge have 
been presented based on data taken with WFI at the 
ESO/MPG 2.2m telescope and with SOFI at the ESO 
NTT, respectively. This dataset has been supplemented 
by the NICMOS photometry from Paper II, and by pub- 
lic 2MASS data. The near-IR CMD has been statistically 
decontaminated from the disk stars in the foreground, pro- 
ducing a very clean CMD of the bulge. To ensure great 
statistical significance also for the upper RGB, the SOFI 
data have been combined with 2MASS data on a larger 
field including the SOFI field itself, and the 2MASS CMD 
has also been statistically decontaminated from the fore- 
ground stars. 

By combining the best near-IR data in each corre- 
sponding luminosity range (2MASS, SOFI, and NICMOS) 
a luminosity function of the Galactic bulge has been con- 
structed that extends over ~ 15 magnitudes, from a few 
bright AGB stars above the RGB tip, down to almost the 
bottom of the main sequence, corresponding to stars of 
~ 0.15 Mq. This is the most extended and complete lu- 
minosity function so far obtained for the galactic bulge, 
hence for galactic spheroids in general. 

Combining WFI optical data and 2MASS near-IR 
data, a disk decontaminated (K, V — K) CMD was con- 
structed that includes 503 stars brighter than Mr = —4.5. 
Together with empirical RGBs of globular clusters of 
known metallicity used as templates, this fairly large sam- 
ple of bulge stars is used to determine the bulge metallicity 
distribution. It is found to peak at [M/H]~ —0.1, with a 
fairly sharp edge just above solar metallicity, and a low- 
metallicity tail that does not appreciably extend below 
[M/H]~ —1, quite similar to the spectroscopic result of 
McWilliam & Rich (1994), but relative to a much larger 
sample. Such a distribution contains somewhat less metal 
poor stars than predicted by the closed box model of chem- 
ical evolution (e.g., Rich 1990), and indicates that the clas- 
sical G-dwarf problem may affect also the Galactic bulge, 
as known since a long time for the disk (e.g., Tinsley 1980). 
Integrated light studies indicate that a G-dwarf problem is 
also present in early type galaxies, as a closed-box model 
metallicity distribution will result in too low values of the 
Mg2 index (Greggio 1997). Yet, the reasons for both disks 
and spheroids sharing the same problem are likely to be at 
least in part quite different. For the Galactic disk (and per- 
haps disks in general) pre-enrichment by the older bulge is 
likely to set a floor metallicity at [Fe/H]~ —1.0 (Rcnzini 



2002). Clearly, this kind of pre-enrichment cannot be in- 
voked for the bulge itself. 

In the closed-box model one assumes the whole bary- 
onic mass to be already assembled at t = 0, while by 
adopting a gradual built up the G-dwarf problem can be 
eliminated. In the case of the Galactic disk this solution of 
the problem has taken the name of infall (e.g. Pagel 1997, 
and references therein). Since the star- forming ISM (ei- 
ther in a disk or in a proto-spheroid) is likely to be built 
up gradually, then the closed-box model would not ap- 
ply. This situation is indeed established both in the case 
of continuous gas infall forming disks, as well as in the 
case of early, merging-driven starbursts building galactic 
spheroids. In turn, starbursts are likely to drive metal- 
rich galactic winds, thus discontinuing star formation be- 
fore the gas is fully turned into stars. Therefore, the bulge 
likely departed from the closed-box approximation both 
for having being built gradually, albeit rapidly, and for 
having ejected gas and metals near the climax of its star 
formation activity. 

The near-IR CMD of the Galactic bulge, carefully de- 
contaminated from foreground stars, has been used to 
show that the bulge is virtually coeval with the Galactic 
halo. The region in the CMD above the main sequence 
turnoff of the bulge is so devoid of stars that no trace 
can be found of a population of bulge stars significantly 
younger than the main old component. This indicates that 
no appreciable new stars formed in the bulge, or were 
added to it, following the intense starburst activity that 
turned ~ 1O 1O M0 of gas into the bulge stars, some 10 Gyr 
ago. 

Using both optical and near-IR data, the RGB tip of 
the bulge is clearly identified, and its systematics in the 
various bands is established. Contrary to the behavior at 
lower metallicities, in the metallicity range spanned by 
the bulge stars (— l<[Fe/H]<0) the /-band luminosity of 
the RGB tip is not constant, but decreases with increasing 
metallicity. Conversely, the RGB tip luminosity is found to 
be constant in the H band, while it increases with metal- 
licity in the K band. This behavior is naturally accounted 
for by a combination of increasing blanketing which is 
stronger at short wavelengths (i.e. in the / and J band), 
with the bolometric luminosity of the RGB tip which in- 
creases with increasing metallicity. The metallicity distri- 
bution of galactic spheroids is therefore expected to affect 
the use of the surface brightness fluctuation method of 
distance determinations (Tonry & Schneider 1988). As a 
final remark, we would like to note that a group of stars 
is found at luminosities exceeding the RGB tip, and their 
number is roughly consistent with the expectation scal- 
ing from old, metal rich globular clusters. Their number 
is consistent with the purely old age for the bulge popu- 
lation as derived from the turnoff region. 

Worth mentioning are also some limitations of the 
present study. The location of the peak of the metallicity 
distribution coincides with the metallicity of the clusters 
NGC 6528 and NGC 6553, but in turn the metallicity it- 
self of these clusters remain quite uncertain. Systematic, 
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homogeneous, high spectral resolution studies of stars in 
these clusters and in the field are needed to properly match 
the cluster and field star metallicity scales. When such 
data will become available, it will be easy to re-derive the 
metallicity distribution of the bulge. Finally, we regret to 
have been unable to find the origin of the larger HB and 
RGB spread exhibited by the bulge CMD, compared to 
simulations that are meant to include all intrinsic spreads 
of the real bulge population (distance, metallicity, and red- 
dening) . 

In closing, we would like to note that the Milky Way 
galaxy is a rather uncospicuous late-type spiral, member 
of a poor group of galaxies located quite away from major 
density peaks in the distribution of galaxies. Yet, virtually 
its whole spheroidal component, from the outer metal- 
poor halo all the way to the central metal-rich bulge, is 
~ 1 Hubble time old. There is now conclusive evidence 
that the stellar populations of galactic spheroids in gen- 
eral, elliptical and bulges, are also ~ 1 Hubble time old 
(e.g. Renzini 1999, for an extensive review), while most 
such spheroids appear to be well in place and fully assem- 
bled already at z ~ 1, having passively evolved since at 
least z ~ 2.5 (Cimatti et al. 2002). If the stellar mass in 
galactic spheroids amounts to 50-70% of the total stellar 
mass in the local universe (Fukugita et al. 1998), all this 
evidence indicate that a major fraction of all stars formed 
at high redshift, in deep potential wells bound to become 
the galactic spheroids of today. If such build up of galaxies 
was promoted by the hierarchical merging of dark matter 
halos, it took place at a time and at a rate that no practi- 
cal rendition of the hierarchical merging paradigm has yet 
succeded to reproduce (Cimatti et al. 2002 and references 
therein). 
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